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PREFACE 


This  report  is  published  to  provide  coastal  engineers  with  a  formulation 
and  a  set  of  nomographs  for  determining  the  breaking  wave  characteristics, 
such  as  breaking  wave  height,  depth  of  breaking,  and  angle  of  breaking  wave 
with  a  straight  shoreline,  as  functions  of  the  deepwater  wave  characteristics: 
wave  height,  wave  period,  and  wave  angle.  This  formulation  is  necessary  to 
determine  the  littoral  drift  transport;  however,  to  obtain  such  results,  a 
review  of  nonlinear  wave  transformation  is  presented.  A  "hybrid"  wave  ap¬ 
proach  based  on  linear  (or  Stokes  third  order)  and  cnoidal  waves  is  proposed 
as  the  best  theory  from  available  experimental  data.  The  work  was  carried 
out  under  the  coastal  structures  program  of  the  Coastal  Engineering  Research 
Center  (CERC) . 

This  report  was  prepared  by  Bernard  Le  Mehaute,  Professor  and  Chairman, 
and  John  D.  Wang,  Associate  Professor,  Ocean  Engineering,  Rosenstiel  School 
of  Marine  and  Atmospheric  Science,  University  of  Miami,  Miami,  Florida, 
under  CERC  Contract  No.  DACW72-79-C-0005 . 
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gation  and  for  the  review  of  the  final  report. 
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CONVERSION  FACTORS,  U.S.  CUSTOMARY  TO  METRIC  (SI)  UNITS  OF  MEASUREMENT 

U.S.  customary  units  of  measurement  used  in  this  report  can  be  converted  to 
metric  (SI)  units  as  follows: 


Multiply 

by 

To  obtain 

inches 

25.4 

millimeters 

2.54 

centimeters 

square  inches 

6.452 

square  centimeters 

cubic  inches 

16.39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

square  yards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  meters 

miles 

1.6093 

kilometers 

square  miles 

259.0 

hectares 

knots 

1.852 

kilometers  per  hour 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

millibars 

1.0197  x  10"3 

kilograms  per  square  centimeter 

ounces 

28.35 

grams 

pounds 

453.6 

grams 

0.4536 

kilograms 

ton,  long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angle) 

0.01/45 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins1 

*To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  readings,  use 


formula:  C  =  (5/9)  (F  -52). 

To  obtain  Kelvin  (K)  readings,  use  formula:  K  =  (5/9)  (F  -32)  +  273.15. 
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SYMBOLS  AND  DEFINITIONS 
A  elliptic  function 

B  elliptic  function 

C  wave  phase  speed 

E  complete  elliptic  integral  of  the  second  kind 

H  wave  height 

unrefracted  deepwater  wave  height 
K  complete  elliptic  integral  of  the  first  kind 

Kr  refraction  coefficient 

Ks  shoaling  coefficient 

L  wavelength 

la  cnoidal  wavelength  parameter 

Lj  linear  wavelength 

Q  littoral  transport  rate 

S  beach  slope 

T  wave  period 

U  Ursell  parameter  *  HL2/d3 

V(u,v)  particle  velocity 
d  Stillwater  depth 

ffj  cnoidal  shoaling  function 

g  gravity 

k  wave  number 

p  pressure 


;-y  cartesian  ordinates 

i  angle  of  wave  incidence 

i  free-surface  elevation 

density 

i(x,z,t)  velocity  potential 

>,o  subscripts  which  refer  to  breaking  or  deepwater  values,  respectively 
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TRANSFORMATION  OF  MONOCHROMATIC  WAVES  FROM  DEEP 
TO  SHALLOW  WATER 

by 

Bernard  Le  Mehaute  and  John  D.  Wang 

I.  INTRODUCTION 


An  understanding  of  many  nearshore  phenomena  relies  on  the  ability  to 
predict  the  local  wave  climatology,  given  a  deepwater  wave  description. 

For  example,  a  quantitative  description  of  longshore  sediment  transport  is 
based  on  a  knowledge  of  the  wave  characteristics  in  the  surf  zone.  This 
report  presents  methods  for  determining  the  changes  in  the  characteristics 
of  a  wave  traveling  over  a  variable  bottom  from  deep  water  to  shallow 
water. 

The  acute  sensitivity  of  the  rate  of  littoral  transport  to  wave 
breaking  characteristics  implies  an  accurate  determination  of  these 
characteristics.  The  problem  has  numerous  facets: 

(a)  Given  a  deepwater  unidirectional  monochromatic  wave,  what 
are  the  breaking  wave  angle,  depth  of  breaking,  breaking  wave  height, 
and  related  quantities? 

(b)  -Given  a  multidirectional  deepwater  incident  wave  spectrum, 
what  is  the  distribution  of  breaking  wave  characteristics  and  the 
"equivalent”  monochromatic  wave  used  to  determine  the  littoral  drift? 

(c)  How  should  a  synoptic  wave  climatology  be  treated  in  order 
to  determine  the  rate  of  littoral  drift  and  related  quantities? 

Only  the  first  problem  is  addressed  in  this  report.  The  relevant 
literature  is  reviewed,  and  a  new  hybrid  wave  theory  is  proposed  to 
determine  wave  breaking  characteristics  on  a  sloped  plane  beach. 

II,  NONLINEAR  WAVE  TRANSFORMATION 
1 .  Nonlinear  Wave  Shoaling- 

It  is  generally  assumed  that  the  wave  motion  over  a  gentle  slope 
is  the  same  as  that  on  a  horizontal  bottom,  and  that  there  is  no  re¬ 
flection  nor  wave  profile  deformation.  The  wave  motion  is  then 
determined  so  that  the  rate  of  transmission  of  energy  or  energy  flux 
is  constant  over  varying  depth. 

The  average  energy  flux  through  a  vertical  plane  of  unit  width 
perpendicular  to  the  wave  propagation  is 

Fav  -f t  (sz  *  h  I y2)  ud2dt  (,) 
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where 


p  =  density 
t  =  time 
T  =  wave  period 
d  =  water  depth 
n  =  free-surface  elevation 
g  =  gravity  acceleration 
p  =  pressure 
V(u,v)  =  particle  velocity 
z  =  vertical  ordinate 


In  the  general  case,  linear  or  nonlinear,  where  the  flow  motion 
can  be  expressed  by  a  potential  function  <j>(x,z,t),  the  Bernoulli 
equation  yields 

-4’ t  =  sz  +  +  \  y2  (2) 

and  u  =  4>x  so  that  the  energy  flux  becomes 


F 


av 


4)  i,  dxdt 
t  x 


(3) 


in  which  case  <p  can  be  expressed  at  any  order  of  approximation,  such  as 
given  by  a  Stokesian  power  series.  Even  though  classical  solutions  for 
cnoidal  waves  are  irrotational ,  the  potential  function  is  not  expressed 
but  rather  the  solution  for  (n,u,v)  is  given;  therefore,  the  energy  flux 
for  cnoidal  wave  is  determined  from  equation  (1)  where  (V2  =  u2  +  w2) . 

The  results  of  all  the  calculations  pertinent  to  linear  wave 
theory  and  linear  wave  shoaling  are  given  in  Le  Mehaute  ( 1 9 7 ( ) , 

Instead  of  expressing  <j>  at  a  first  order  of  approximation  as  in 
the  linear  wave  theory,  <t>  is  expressed  at  a  higher  order  in  equat  ion 
(3),  the  shoaling  coefficient  Kg  =  H/H0  becomes  not  only  a  function 
of  d/L  or  d/L0  but  also  a  function  of  the  deepwater  wave  steepness, 

H0/L0 • 


r 


This  calculation  has  been  performed  at  a  third  order  of  approxi¬ 
mation  (Le  Mehaute  and  Webb,  1964),  and  the  fifth  order  of  approximation 
(Koh  and  Le  Mehaute,  1966)  based  on  the  third-order  solution  and  fifth- 
order  solution  for  a  Stokesian  wave  as  developed  by  Skjelbreia  and 
Hendrickson  (1960).  The  first  definition  of  Stokes  for  the  phase  velocity 
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is  used;  i.e.,  the  average  horizontal  water  particle  velocity  over  a 
wavelength  is  zero.  The  results  of  such  investigation  are  presented 
in  Figures  1  and  2. 

The  correction  AH  due  to  nonlinear  effects  never  exceeds  5 
percent  and  is  more  commonly  of  the  order  of  1  percent.  These  in¬ 
vestigations  show  that: 

(a)  The  nonlinear  shoaling  coefficient  is  initially  less  than 
the  linear  coefficient  when  d(  L0  5  0.4,  then  becomes  larger  toward 
shallow  water  until  the  wave  Dreaks. 

(b)  The  Stokesian  power  series  is  not  uniformly  convergent,  i.e., 
the  function  of  d/L  of  higher  order  tends  toward  infinity  when  c/L  tends 
to  small  values.  Therefore,  the  "best"  order  of  approximation  is  not 
necessarily  the  highest  order.  For  relatively  deep  water  d/L  >  0.25,  the 
fifth  order  of  approximation  would  be  the  best  insofar  as  wave  height 
transformation  is  concerned;  for  very  shallow  water  d/L  <  0.01,  the 
linear  theory  would  be  best.  In  the  intermediate  range  the  third-order 
theory  would  be  best,  and  therefore  should  be  preferred  overall  because 
of  its  range  of  applicability. 

The  second  definition  of  Stokes  for  the  phase  velocity  can  also  be 
used;  the  average  momentum  over  a  wavelength  is  zero  by  addition  of  a 
uniform  motion.  Yamaguchi  and  Tsuchiya  (1976)  indicate  that  the  results 
yield  slightly  larger  values,  at  most  a  7-percent  increase  for  the 
shoaling  coefficient,  than  the  results  obtained  by  Le  Mehaute  and  Webb 
(1964). 

The  principle  of  conservation  of  energy  flux  has  also  been  applied 
to  a  cnoidal  wave,  and  like  the  Stokesian  wave  the  resuits  depend  on  the 
order  of  approximation  and  the  definition  of  phase  velocity.  All  these 
investigations  on  cnoidal  waves  are  based  on  ai.  energy  flux  such  as  ex¬ 
pressed  by  equation  (1).  Masch  (1964)  was  the  f •  •  s t  to  d»al  with  this 
subject;  however,  his  wave  theory  is  not  consistent,  even  erroneous, 

(in  the  table  of  functions  used  by  Masch  in  the  shoaling  of  cntwdal 
wave,  the  water  depth  below  MWL  should  be  substituted  by  ht,  the  water 
depth  under  trough) ,  and  the  results  are  presented  in  a  form  which  is 
difficult  to  use.  The  relation  to  deepwater  wave  and  sinusoidal  theory 
is  not  discussed  and  no  attempt  is  made  to  follow  the  shoaling  of  a 
specific  wave. 

A  significant  contribution  to  the  shoaling  of  cnoidal  waves  is 
given  by  Iwagaki  (1968).  Iwagaki  treats  the  case  of  an  approximate 
solution  of  cnoidal  wave  in  which  he  used  the  second  definition  of 
phase  velocity  as  given  by  Laitone  (1961).  The  approximation  is  on  the 
value  of  the  elliptic  integral  which  is  replaced  by  a  simple  function 
of  empirical  coefficients.  Iwagaki  shows  that  this  simplification 
actually  covers  a  wide  range  of  cases  and  allows  him  to  simply  in¬ 
vestigate  the  shoaling  of  what  he  calls  "hyperbolic  waves."  When  the 
energy  flux  in  deep  water  (as  computed  using  small -amplitude  theory) 
is  equated  to  the  energy  flux  in  shallow  water,  described  by  first- 
order  hyperbolic  waves,  Iwagaki  obtains 
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According  to  Iwagaki  (1968),  this  theory  yields  sufficiently  accurate 
results  for  Ursell  parameter  U  >  47.  However,  as  pointed  out  by 
Svendsen  (1974),  the  theory  of  Iwagaki  deserves  to  be  regarded  as  a 
practical  solution  to  second-order  cnoidal  waves  when  the  deepwater 
wave  steepness  is  smaller  than  0.02  and  the  relative  water  depth  is 
smaller  than  0.05.  The  matching  of  the  Iwagaki  hyperbolic  wave  with 
the  third-order  Stokesian  wave  is  shown  in  Figure  3. 

The  shoaling  of  the  true  cnoidal  wave  has  been  investigated  by 
Svendsen  and  Brink-Kjaer  (1972),  Svendsen  (1974),  and  Svendsen  and 
Hansen  (1977).  They  also  give  H/Hn  as  function  of  d/L0  and  H0/L0 
(Fig.  4)  and  a  computer-printed  table.  It  can  then  be  shown  that  for 
large  values  of  Ursell  parameters  the  shoaling  coefficient  Ks  -*•  d~l 
instead  of  d“V4  as  given  by  the  Green  law  (long  wave  linear  theory) . 
Concurrently,  Shuto  (1974)  arrives  at  very  similar  results. 

Yamaguchi  and  Tsuchiya  (1976)  also  carry  out  the  same  calculation 
based  on  the  two  definitions  of  the  Stokes  wave  velocity  for  the  cnoidal 
theory  of  Laitone  (1961)  and  that  of  Chappelear  (1962).  However,  an 
arithmetic  error  has  been  found  in  the  laitone  theory  (Le  Mehaute,  1968) . 

2.  Comparison  and  Matching  Between  Various  Theories. 

As  a  wave  propagates  from  deep  water  to  shallow  water  it  is 
theoretically  possible  to  determine  the  variation  of  wave  height,  wave¬ 
lengths,  etc.  This  could  be  done  by  applying  the  principle  of  con¬ 
servation  of  energy  flux  to  either  the  linear  wave  or  the  nonlinear 
Stokesian  wave,  or  the  cnoidal  and  solitary  wave.  Since  a  Stokesian 
wave  rather  applies  in  deep  water,  the  transformation  of  water  wave 
should  be  followed  with  that  theory  for  the  largest  value  of  relative 
depth  d/L0  and  then  switched  to  the  cnoidal  theory  when  d/Lo  becomes 
small.  However,  such  a  scheme  implies  that  the  theories  can  be  matched 
continuously,  but  there  is  a  priori  no  reason  why  the  ratio  H/H^  should 
be  the  same  for  the  value  d/L0  which  corresponds  to  the  limit' of  validity 
of  both  theories.  On  the  other  hand,  if  the  wave  heights  are  matched, 
then  the  energy  flux  will  present  a  discontinuity  (Fig.  5).  The  signifi¬ 
cant  feature  is  that  the  cnoidal  wave  height  grows  faster  with  decreasing 
depth,  although  at  intermediate  depth  its  value  is  up  to  10  percent  less 
than  predicted  by  a  Stokesian  theory.  Waves  with  wave  steepness  larger 
than  2  to  3  percent  will  break  at  a  depth  where  the  cnoidal  wave  height 
is  only  slightly  larger  than  that  of  a  Stokesian  wave.  Waves  with  small 
wave  steepness,  however,  such  as  swells,  reach  much  smaller  depth  before 
they  break  and  consequently  a  major  part  of  their  shoaling  process  is 
governed  by  the  cnoidal  wave  theory.  For  these  waves,  the  two  theories 
such  as  the  Stokesian  (first  order  or  linear  theory)  and  cnoidal  wave 
at  a  second  order  will  yield  significantly  different  results. 


H/M 
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d/L* 

Figure  5.  Matching  Stokesian  (first  order)  and  cnoidal 
wave  theory  (Svendsen  and  Hansen,  1977) . 


These  results  (Fig.  5)  show  that  no  continuous  transition  is  pos¬ 
sible  between  the  two  theories.  This  means  that  it  is  not  possible  to  find 
a  value  of  the  water  depth,  d,  where  the  curves  for  the  two  theories  fit 
smoothly  together.  If  the  Stokesian  theory  is  used  in  deeper  water  and 
changed  to  a  cnoidal  theory  when  the  wave  enters  shallow  water,  there 
will  be  a  discontinuity  in  the  variation  of  either  wave  height  or  wave¬ 
length,  or  both,  depending  on  which  water  depth  is  chosen  for  the  switch. 

Of  course,  the  same  will  appear  for  all  other  quantities  such  as  particle 
velocities,  pressure,  etc.,  and  the  rate  of  change  of  these.  Svendsen 
(1974)  shows  that  the  limit  of  applicability  of  the  cnoidal  theory  is 
d/L0  <  0.1193  when  H  is  small.  Koh  and  Le  Mehaute  (1966)  also  showed 
that  the  limit  of  applicability  of  the  fifth-order  Stokesian  wave  theory 
is  d/LQ  >  0.10  when  II/LQ  =  0.05  and  d/L0  >  0.13  when  H/L0  =  0.10  (see 

Fig.  2). 

There  is  a  large  difference  between  Stokesian  and  cnoidal  wave 
between  d/L0  equal  0.1  and  0.3.  In  this  region  no  known  wave  theory 
fits  very  well.  It  could  have  been  expected  that  a  higher  order 
Stokesian  theory  would  be  the  answer,  but  the  investigation  by  Koh 
and  Le  Mehaute  (1966)  shows  that  when  d/L0  decreases  the  fifth-order 


approximation  represents  an  even  worse  approximation  than  the  third 
order.  Similarly,  it  is  found  that  second-order  cnoidal  theory  is 
worse  than  first -order  cnoidal  theory  for  large  wave  steepness.  This 
is  inherent  to  the  point  that  both  cnoidal  and  Stokesian  power  series 
expansion  in  terms  of  the  small  parameters  h/d  and  H/L  respectively 
are  nonuniformly  converging  series  since  the  functions  of  d/L  attached 
to  each  power  term  blow  up  when  d/L  tends  toward  small  values. 


It  is  interesting  that  Yamaguchi  and  Tsuchiya  (1976)  found  that 
the  shoaling  coefficient  given  by  Le  Mfehaute  and  Webb  (1964)  (first 
definition  of  Stoke' s  phase  velocity)  almost  coincides  with  the 
shoaling  coefficient  obtained  from  cnoidal  theory  developed  by 
Chappelear  (1962)  (second  definition). 


Shuto  (1974)  attempted  to  make  a  synthesis  of  all  these  theories 
in  a  simple  and  practical  form  by  empirically  matching  these  solutions. 
Subsequently,  he  proposes  the  following  law  for  practical  purposes: 
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These  equations  seem  to  be  the  most  realistic  to  remember  from  all 
the  theoretical  approaches.  In  the  range  where  both  cnoidal  and  third- 
order  Stokesian  theory  apply,  the  values  of  the  shoaling  coefficient 
are  very  close  to  each  other  as  shown  in  Figure  6  (Flick,  1978). 


1 


d/Lo 


Figure  6.  Comparison  between  Stokesian  third  order  and 
cnoidal  shoaling  coefficient  with  experiments 
(from  Flick,  1978) . 
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Interestingly,  the  use  of  the  linear  wave  theory  to  evaluate  the 
value  of  the  shoaling  coefficient  extends  much  beyond  the  formal 
validity  of  this  infinitesimal  wave  theory.  Similarly,  the  value  of 
the  shoaling  coefficient  given  by  the  Stokesian  wave  theory  extends 
into  the  area  where  the  cnoidal  theory  fits  best.  This  is  due  to  the 
fact  that,  the  shoaling  coefficient  being  the  ratio  of  wave  height 
H/H0  only,  the  increase  in  free-surface  elevation  under  the  crest  is 
partly  balanced  by  the  increase  of  free-surface  elevation  under  the 
wave  trough.  However,  that  the  linear  wave  theory  applies  for  the 
shoaling  coefficient  does  not  mean  that  all  wave  characteristics 
(wavelength,  velocity  components,  pressure,  acceleration)  follow  the 
same  principle;  after  the  local  wave  height  is  obtained,  all  other 
wave  characteristics  are  determined  by  the  appropriate  theory. 

3.  Comparison  Between  Theory  and  Experiment • 

A  relatively  large  number  of  experiments  have  attempted  to 
verify  shoaling  laws;  all  have  been  conducted  in  laboratory  wave  flumes 
with  waves  generated  by  wave  paddle.  Most  of  these  experiments  suffer 
lack  of  accuracy  because  they  were  either  done  at  too  small  a  scale 
and  were  subsequently  subjected  to  significant  scale  effects  such  as 
large  viscous  damping  experiments  (Iversen,  1951),  or  the  wave  paddle 
generated  not  only  monochromatic  waves  but  harmonic  components 
(solitons)  which  introduced  significant  error  and  scattering 
(Eagleson,  1956;  Iwagaki,  1968). 

There  is  actually  considerable  controversy  whether  waves  of 
steady-state  profile  exist,  as  demonstrated  by  Dubreuil-Jacotin  (1934). 

Theorists  Benjamin  and  Feir  (1967)  and  experimentalist  Galvin  (1970) 
postulate  that  the  disintegration  of  finite  amplitude  monochromatic 
wave  occurs  in  deep  water  even  on  horizontal  bottom.  There  are  as  many 
theoreticians  who  assume  that  a  steady-state  profile  does  exist  as  there 
are  experimentalists  who  do  not  notice  the  "creations"  of  solitons. 

Use  of  a  formulation  developed  by  Mei  and  Le  Mehaute  (1966), 

Peregrine  (1967),  and  Madsen  and  Mei  (1969)  indicates  that  for  a 

sufficiently  abrupt  change  in  water  depth,  both  a  solitary  wave  and  a 

cnoidal  wave  disintegrate  into  multiple  crests.  These  results  have 

been  obtained  numerically  and  verified  experimentally.  However,  over 

a  relatively  gentle  beach,  the  wave  period  remains  constant  between 

deep  water  and  shallow  water  and  no  disintegration  takes  place.  Dis-  j 

integration  takes  place  when  the  wave  arrives  on  a  reef.  It  seems  i 

natural  to  assume  that  the  difference  between  these  two  observations  j 

is  due  to  the  difference  in  bottom  slope.  Benjamin  and  Feir  (1967)  j 

show  that  waves  are  unstable  if  kd  >  1.4;  however,  experiments  by 

Flick  (1978)  indicate  that  kd  can  be  much  larger  without  evidence  of 

wave  disintegration  or  spectral  smearing. 

It  is  commonly  accepted  that  a  monochromatic  wave  arriving  on  a 
rapid  change  of  depth  (in  diffraction  zone)  gives  rise  to  at  least  a 
doubling  of  crests.  Such  phenomenon  is  due  to  the  nonlinear  con¬ 
vective  effects.  Iwagaki  and  Sagai  (1971)  also  investigated  the 
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deformation  of  long  waves  over  a  gentle  slope  using  the  nonlinear  long 
wave  theory  and  power  series  expansions.  They  found  the  shoaling 
coefficient  to  be  a  function  of  beach  slope  when  S  >  .01.  The  steeper 
the  slope  the  smaller  the  shoaling  coefficient,  a  fact  which  can  be 
attributed  to  partial  reflection.  In  fact,  due  to  friction  effect,  the 
ratio  H/H0  for  a  given  value  of  d/L0  decreases  instead  of  increases 
(Sawaragi,  Iwata,  and  Masayashi,  1976)  . 

The  first  reliable  experiments  were  conducted  by  Brink-Kjaer  and 
Jonsson  (1973)  and  Flick  (1978) .  Figures  7  and  8  show  results  for  different 
values  of  H0/LQ.  Flick  separates  the  first,  second,  and  third  harmonics 
from  his  wave  data  and  is  subsequently  able  to  give  a  reliable  ex¬ 
perimental  shoaling  coefficient.  Flick  compares  his  results  with 
Le  Mehaute  and  Webb  (1964)  (third-order  Stokesian)  and  also  with  the 
cnoidal  solution  of  Svendsen  and  Brink-Kjaer  (1972)  in  shallow  water 
(see  Fig.  6) . 

The  shoaling  coefficient  of  a  hyperbolic  wave  is  also  fairly  well 
verified  by  Iwagaki  (1968)  who  gives  results  very  close  to  the  two 
mentioned  above. 

Svendsen  and  Hansen  (1977)  compared  the  shoaling  of  cnoidal  wave 
with  a  set  of  careful  experiments  and  claimed  that  other  experimenters 
(Wiegel,  1950,  Iversen,  1951;  Eagleson,  1956)  carried  out  their  ex¬ 
periments  on  too  steep  a  slope  for  the  shoaling  theory  to  be  valid. 
Furthermore,  they  calculate  the  damping  due  to  viscous  friction, 
obviously  important  on  a  gentle  slope.  Svendsen  and  Hansen  concluded 
that  if  the  wave  height  at  depth  d/L0  =  0.10  is  matched  between  cnoidal 
and  linear,  rather  than  the  energy  flux,  the  cnoidal  theory  predicts 
the  shoaling  quite  well,  even  close  to  breaking  with  small  deepwater 
wave  steepness  H0/L0  <  3  to  4  percent  but  not  beyond.  Consistently, 
with  all  theories,  the  wave  just  before  breaking  suddenly  peaks  up  very 
rapidly  (Le  Mehaute,  1971).  In  this  range  of  values,  all  shoaling 
theories  (third  Stokes,  cnoidal  and  hyperbolic)  tend  to  slightly  under¬ 
estimate  the  value  of  the  shoaling  coefficient.  Subsequently,  the  cal¬ 
culated  breaking  wave  height  tends  to  be  underestimated.  The  linear 
wave  theory  underestimates  the  breaking  wave  height  most  significantly, 
sometimes  by  a  factor  of  almost  2  (Fig.  9). 

It  is  pertinent  to  remember  that  (a)  the  shoaling  coefficient 

given  by  the  linear  theory  is  valid  beyond  the  limit  generally  con¬ 

sidered  as  valid  for  a  linear  theory,  and  (b)  the  shoaling  coefficient 
given  by  third-order  Stokesian  wave  is  fairly  well  verified  ex¬ 
perimentally  and  actually  very  ’ose  to  the  value  given  for  the  cnoidal 
wave,  even  though,  as  in  the  c ?  .e  of  the  linear  wave,  free-surface 
profile,  pressure,  velocity,  and  acceleration  could  be  significantly 
different . 

In  general,  the  linear  theory  can  be  applied  throughout  from  deep 

water  to  shallow  water  and  then  the  linear  breaking  wave  height  is  multi¬ 

plied  by  a  coefficient  function  of  the  beach  slope  (Koh  and  Le  Mehaute, 
1966).  After  the  wave  height,  H,  is  determined  as  a  function  of  the 
deepwater  wave  height,  HQ,  and  wave  period,  T  (or  deepwater  wavelength 
L0),  all  other  shallow-water  characteristics  f free-surface  profile, 
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particle  velocity,  acceleration,  and  pressure)  follow  by  application 
of  one  of  the  classical  wave  theories  within  the  accuracy  which  is 
determined  for  the  chosen  theory. 


4.  Nonlinear  Wave  Refraction. 

It  has  been  shown  previously  how  to  determine  the  shoaling  co¬ 
efficient,  Ks  =  H/H0,  when  the  wave  arrives  perpendicular  to  the 
bottom  contour.  This  discussion  deals  with  the  refraction  coefficient, 
Kr.  Refraction  occurs  when  a  wave  arrives  at  an  angle  a  with  bottom 
contours;  then  H/HQ  =  Ks  Kr.  For  a  straight  parallel  contour  Snell's 
law  becomes 

—  =  constant  (6) 

sin  «  ' 

which  applies  whether  the  wavelength,  L,  is  expressed  by  a  linear  theory 
or  not. 


21 


Esssa 


•im/wvsw&mr. 


which  also  applies  for  nonlinear  as  well  as  for  linear  theory.  The 
subscript  o  refers  to  deepwater  wave  characteristics. 

In  many  cases,  the  refraction  method  provides  a  reasonably 
accurate  measure  of  the  changes  waves  undergo  on  approaching  a  coast. 

However,  if  the  angle  of  a  wave  ray  with  the  bottom  contour  is  large 
(i.e.,  larger  than  70°),  minor  error  in  the  value  of  the  incident 

angle  leads  to  a  large  error  in  direction  angle  a  in  shallow  water.  ! 

Also,  accuracy  as  far  as  height  changes  are  concerned  cannot  be  ex-  j 

pected  where  bottom  slopes  are  steeper  than  1/10.  No  strict  limit  .* 

has  been  set,  but  the  accuracy  of  wave  heights  derived  from  orthogonals  \ 

that  bend  sharply  is  questionable.  In  short,  refraction  coefficients  ^ 

which  are  quite  different  from  unity,  such  as  Kr  <  0.5  and  Kr  >  1.5,  j 

must  be  doubted  (Whalin,  1971).  •  » 

Nonlinear  effects,  having  an  effect  on  wavelength,  phase  and  j 

group  velocity  and  energy  flux,  subsequently  have  an  effect  on  wave  *  . 

refraction.  This  problem  has  been  examined  by  Chu  (1S75)  who  used  a  I 

mix  of  three  theories,  i.e.,  the  first-order  cnoidal  theory  of  Korteweg  j 

and  DeVries  (1895),  the  second-order  hyperbolic  wave  of  Iwagaki  (1968),  l 

and  the  Stokes  third-order  wave  as  given  by  Le  Mehaute  and  Webb  (1964) ,  J 

which  led  to  some  inconsistencies  in  approximations.  Skovgaard  and  ' 

Petersen  (1977)  used  instead  the  first-order  cnoidal  theory  of  Svendsen  [ 

(1974)  and  the  stream  function  wave  theories  of  Dean  (1970). 

Theoretically,  it  is  possible  to  express  phase  velocities  as  a  ' 

function  of  the  relative  wave  heights  from  nonlinear  wave  theories. 

For  example,  the  deepwater  wavelength  at  a  third  order  Stokesian 
approximation  and  the  breaking  wavelength  by  a  cnoidal  or  hyperbolic 
wave  theory  can  be  conveniently  expressed.  However,  it  is  interesting 
that  due  to  deformation  of  wave  profile  on  a  sloped  bottom,  the  simple 
linear  theory  has  been  verified  (experimentally)  quite  well  (Ippen, 

1966) .  Wavelengths  given  by  linear  and  cnoidal  theories  are  compared 
in  Figure  10.  Although  the  cnoidal  theory  predicts  wave  height  well  up 
to  breaking,  it  overpredicts  wavelengths  significantly.  Cnoidal  theory, 
in  fact,  predicts  an  increase  in  wavelength  for  a  decrease  in  depth 
when  the  relative  height,  H/d,  is  sufficiently  large.  This  increase 
is  not  reflected  by  known  data  (Ippen,  1966)  which  are  fitted  quite 
well  by  linear  theory  (Fig.  11). 

III.  BREAKING  WAVE  CHARACTERISTICS  ON  A  SLOPED 
PLANE  BEACH 

1 .  Review  of  Previous  Work  . 

The  determination  of  longshore  currents  and  sediment  transport 
depends  crucially  on  the  characteristics  of  the  breaking  wave  field. 

The  wave  energy  transport,  or  energy  flux,  is  of  particular  importance 
such  that  accurate  determination  of  wave  height,  wavelength,  depth  at 
breaking,  and  breaking  wave  angle  becomes  essential. 

! 

This  section  deals  with  the  practical  aspects  of  determining  [ 

the  breaking  wave  characteristics  when  certain  deepwater  character-  j 

istics  are  given.  The  objective  is  to  derive  and  present  results  [ 
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Figure  10.  Wavelength  as  a  function  of  depth  pre¬ 
dicted  by  different  wave  theories. 


d/L 

Figure  11.  Transformation  of  wavelengths  (.from  Ippen, 
1966;  used  with  the  permission  of  McGraw- 
Hill  Book  Company). 
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consistent  with  present  knowledge  and  in  a  readily  usable  form. 

The  general  problem  would  require  the  determination  of  the 
shoaling  and  refraction  of  a  multidirectional  wave  spectrum  from  deep 
water  over  a  randomly  varying  bottom  topography  until  breaking  occurs. 
Although  such  an  analysis  is  possible,  it  is  much  too  complicated  and 
would  have  to  be  dealt  with  on  a  case  by  case  basis  using  either 
manual  or  computer  methods. 


A  significant  and  useful  simplification  is  achieved  by  assuming 
the  bottom  to  be  a  uniformly  sloped  plane.  This  allows  bottom  vari¬ 
ations  to  be  described  by  a  single  parameter,  i.e.,  the  bottom  slope 
S.  The  refraction  process  is  then  described  globally  by  Snell's  law. 

This  discussion  deals  only  with  monochromatic  waves  under  the  usual 
assumption  that  the  wave  period  remains  constant,  and  that  friction 
and  reflection  are  ignored. 

To  obtain  accurate  predictions  it  is  necessary  to  have  a  wave  theory 
which  is  applicable  up  to  the  point  of  breaking.  A  lack  of  knowledge 
of  the  actual  breaking  process  requires  the  use  of  an  empirical 
breaking  criterion  to  determine  the  point  of  breaking.  A  study  by 
Le  Mehaute  and  Koh  (1967)  evaluated  the  Stokes  first-,  third-,  and 
fifth-order  theories  and  compared  the  Miche  (1944)  breaking  criterion 
with  an  empirically  derived  equation.  One  of  these  equations  was 
derived  by  fitting  a  number  of  experimental  data  points  covering  the 
range  0.02  <  S  <  0.2  and  0.002  <  H0/L0  <  0.09.  This  .equation  ex¬ 
plicitly  accounts  for  beach  slope  and  is 


=  0.76  S 


Since  equation  (8)  is  based  on  observed  data  it  takes  into  account 
nonlinear  effects  such  as  wave  height  peak-up  just  before  breaking. 

In  applying  this  equation  to  waves  arriving  at  an  angle  to  the  shore, 

T *  Mehaute  and  Koh  (1967)  corrected  the  bottom  slope  for  the  angle  of 
incidence;  however,  they  neglected  to  replace  the  deepwater  wave 
height  with  its  unrefracted  value. 

Subsequently,  a  new  and  easier  appioach  to  compute  cnoidal 
waves  was  presented  by  Svendsen  (1974),.  Brink-Kjaer  and  Jonsson 
(1973)  showed  that  near  breaking  the  water  depth  is  usually  so  shallow 
that  cnoidal  theory  applies.  Indeed,  it  has  been  found  that  wave 
height  is  described  well  by  cnoidal  theory  in  the  area  close  to  and 
before  breaking. 

In  a  recent  report,  Ostendorf  and  Madsen  (1979)  propose  to  use 
cnoidal  and  linear  Stokes  wave  theories  in  their  respective  areas  of 
applicability.  A  transition  between  the  two  theories  which  assumes 
continuous  variation  of  energy  flux  and  phase  velocities  is  also 
presented.  Ostendorf  and  Madsen  further  suggest  the  use  of  an  em¬ 
pirical  breaking  criterion  which  is  sensitive  to  bottom  slope  and 
depth-varying  wave  parameters,  i.e., 
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(9) 


S  <  0.1 


S  >  0.1 


To  obtain  the  breaking  wave  characteristics,  the  two  offshore 
parameters  (sin  a/C* ,  C4)  must  be  known  where 

a  =  angle  of  incidence 

C  =  wave  phase  speed  C*  =  C/ (gT) 


g  =  gravity  acceleration 


T  =  wave  period 

c4  a  T  ^  {(H°2  n  sin  ct»‘1/4 


d 

n 


=  Stillwater  depth 


2kd  > 
sinh2kdJ 


k  =  2it/l. 


L  =  wavelength 

In  the  deepwater  limit  this  implies  that  wave  height,  H0,  angle 
of  incidence,  a0,  and  wave  period  must  be  known  independently.  The 
solution  requires  an  iteration  process  and  nomographs  are  presented 
to  facilitate  the  operations. 

The  method  for  determining  breaking  wave  characteristics  sug¬ 
gested  by  Ostendorf  and  Madsen  (1979)  has  been  compared  with  experi¬ 
mental  data  (Kamphuis,  1963),  and  it  was  found  that  the  predicted 
breaking  wave  angle  is  too  large,  especially  for  smaller  wave  steep¬ 
nesses  (see  Table  1).  This,  is  easily  explained  when  considering  the 
plot  of  wavelength  transformation  shown  in  Figure  10.  Although 
cnoidal  theory  predicts  wave  bright  well  up  to  breaking,  it  overpre¬ 
dicts  wavelengths  significantly.  Cnoidal  theory,  in  fact,  predicts 
an  increase  in  wavelength  and  therefore,  also  in  wave  angle  when  H/d 
is  sufficiently  large.  This  increase  is,  as  previously  mentioned, 
not  reflected  by  known  data  (Ippen,  1966;  Fig.  11),  which  are  fitted 
quite  well  by  linear  theory.  As  another  consequence,  the  wave  break¬ 
ing  criterion,  again,  a  result  which  is  difficult  to  defend. 

Dean  (1974)  determined  wave  breaking  angles  using  his  stream 
function  theory,  but  with  a  slope- independent  semiempirical  breaking 
criterion.  A  comparison  of  his  results  with  the  experimental  ob¬ 
servations  of  Kamphuis  (1963)  is  also  presented  in  Table  1.  The 
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predictions  are  consistently  too  high,  especially  for  smaller  wave 
steepness,  where  predicted  and  observed  otb  differ  by  a  factor  of 
approximately  2. 


Table  1 .  Comparison  of  measured  and  predicted 
breaking  wave  angles.1 


Ho 

0 

0.0175 

0.04 

0.053 

0.062 

ao 

20° 

o 

o 

Tt- 

o 

60° 

20° 

o 

o 

NJ 

O 

o 

o 

o 

K> 

o 

o 

o 

o 

H 

Kamphuis 

(1963) 

9° 

12° 

8° 

16° 

o 

o 

w—4 

U 

11° 

1 

9 

Ostendorf 
and  Madsen 
(1979) 

18.6° 

24.3° 

12.8° 

24.3° 

13.9° 

26.5° 

14.5° 

% 

Dean 

(1974) 

9.5° 

19° 

22° 

13° 

25° 

14° 

27° 

14.5° 

_ 

o 

00 

CM 

^each  slope,  S  =  .1 


2 .  Solution  Approach. 

To  obtain  reliable  prediction  of  breaking  wave  characteristics, 
this  study  proposes  to  use  cnoidal  theory  to  describe  the  transformation 
of  wave  height  while  wavelength  will  be  transformed  using  either  linear 
wave  theory  or  third-order  Stokes  theory.  The  cnoidal  wavelength  is 
then  considered  as  an  auxiliary  parameter  which  cannot  be  identified  as 
the  physical  wavelength.  Linear  wave  theory  is  simpler  to  use;  however, 
to  retain  some  nonlinear  effect  in  the  transformation  of  wave’ength  the 
third-order  Stokes  theory  is  also  included.  The  wavelength  computed 
using  the  cnoidal  third-order  Stokes  theory  is  shown  in  Figure  11. 

Due  to  the  large  wave  heights  near  breaking,  a  higher  order  approxi¬ 
mation  to  cnoidal  waves  given  by  Iwagaki  (1968)  was  considered.  A  more 
detailed  discussion  of  this  "hyperbolic"  wave  theory  is  given  in  the 
Appendix.  Note  that  this  higher  order  theory  suffers  from  the  same  prob¬ 
lem  of  inhomogeneous  convergence  as,  for  example,  plagues  fifth-order 
Stokes  waves  (Le  Mehaute  and  Koh,  1966),  and  therefore  gives  poorer 
results  than  the  first-order  cnoidal  theory  near  breaking. 

The  computation  and  shoaling  of  cnoidal  waves  have  been  given  by 
Svendsen  (1974)  and  Svendsen  and  Brink-Kjaer  (1972).  It  is  convenient 
to  define  the  parameter 

hl2  ia  ? 

U  =  — y-  =  mK  (m)  (10) 

d 


H 


where 


=  wave  height 

=  cnoidal  wavelength  parameter 


d  =  Stillwater  depth 

K(m)  =  complete  elliptic  integral  of  the  first  kind 


The  dispersion  relationship  may  be  written  as 


r  Vf  - 27r  (1  +  ^ A)  (11) 

0  o 

2  2  2  3F 

The  deepwater  wavelength  L0  =  §7  T  and  A  =  A(m)  =  m  "  *  " 

The  complete  elliptic  integral  of  the  second  kind  is  designated  E. 

Invoking  energy  flux  conservation  between  wave  rays  and  using 
linear  theory  in  deep  water,  the  wave  height  transformation  is  given 
as 

In  equation  (12)  f^  =  f  (U)  =  .  B*2^3 


B  =  B(m)  ^  -  5m  +  2  +  (4m  -  2) 

-  (l  -  m  -  I)2]  (13) 

Equations  (12)  and  (13)  define  the  shoaling  of  cnoidal  waves  and  are 
used  to  determine  the  shoaling  coefficient. 

To  compute  the  wavelength  and  refraction  the  Stokes  wave  theory 
is  used.  Linear  waves  are  described  by  the  dispersion  relationship 

L1  *  u  2*d  n/n 

—  =  tanh  j- —  (14) 

o  L1 

For  the  third-order  Stokes  approximation  from  Le  Mehaute  and  Webb  (1964) 


i  =  I—  tanh 


r.2  "  2tt 


(1  +  Cj  Ap 


where 


>0  a  7 

J  =  27  d  +  ^30) 


l  A3  .  x  =  *Ho 

8  a30  a30  rr- 

30 


(18) 


3  ttH 

3  3  =  L3 

8cosh4kjd  -  8cosh2kjd  +  9 
- - - 

8sinh  kjd 
3(8cosh^k^d  +  1) 
64sinh^k^d 


and 


Refraction  is  described  by  Snell’s  law 


(19) 


(20) 


(21) 


L  _  sing 
Lo  =  sinao 

and  the  refraction  coefficient 


/cosn„ 

(  =  y - EL 

R  T  cosa 


(22) 


(23) 


where  a  =  angle  of  wave  with  shoreline.  Also,  when  refraction  is  included 
the  deepwater  wave  height  in  the  above  expressions  should  be  replaced  by 
its  unrefracted  value. 


Finally,  the  point  of  breaking  is  determined  by  a  modified  form  of 
the  empirical  breaking  criterion  (eq.  8).  To  adapt  the  formula  to  waves 
approaching  at  an  oblique  angle,  the  bottom  slope  and  deepwater  wave 
height  are  replaced  by  S  cosab  and  H0  (cosa0/cosab)^2,  respectively. 

The  breaking  criterion  is  then  obtained  by 


b  „  cl/7  1/7  ("o\ 

-  =  0.76  S  cos  a.  It—) 

o  vo' 


-1/4 
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or 
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Rewriting  equation  (25)  using  H^/Hq  =  KsKr,  where  Ks  is  the  shoaling  co¬ 
efficient  yields 
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K  .  =  0.76 
sb 


iLoS 


o 
4/7 
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-1/4 
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Equation  (?6)  gives  the  form  of  the  breaking  criterion  used  in  this  study. 
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An  explicit  solution  for  the  breaking  wave  characteristics  cannot 
be  obtained  from  equations  (11)  to  (23)  and  equation  (26) .  A  numerical 
solution  is  required  and  it  becomes  important  to  reduce  the  number  of 
independent  parameters  as  much  as  possible.  By  straightforward  manipu¬ 
lation  of  the  equations,  only  three  independent  parameters  need  to  be 
known:  deepwater  wave  steepness,  H0/L0;  beach  slope,  S;  and  deepwater 
incident  wave  angle,  a0. 

A  computer  algorithm  is  constructed  to  solve  the  problem.  The 
basic  process  consists  of  the  following  steps: 

H 

(1)  Give  values  for  j- ,  S,  aQ 

o  i... 


(2)  Assume  Ks=Kr=A=l  and  <  or  l3b| 


(3)  Guess  a  value  for 


(4)  Find 


lit  .  A)  ’*  K  K  %  21,  (l  *  ^*) 

L0  v  Vdb;  S  r  L0  v  \  J 

Hb  Ho  Alb\  _1 

(5a)  Determine  =  Ks  Kr  X  IT  Vl^J 

©-1  Hb 

hb 


^  Ho  fab) 

(Sb)  Determine  IJb  =  Ks  Kr  X 


and  k3b  db  = 


■*© 


-1  H. 


(6a)  Find  j—  =  tanh  kJb  dfe 

10  i  \ 

L3b  UCM 

(6b)  Find  —  =  tanh  k3b  dfe  — - y~ 

3o  \1+X  o3 

(7)  Find  otb  from  equation  (22),  using  wavelength 

ratio  from  step  (6),  and  find  Kr  from  equation 
(16) 
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h,  lib 

(8)  Determine  U  =  x  ^ — 

b  \  o 

(9)  Find  m  and  K(m)  satisfying  equation  (10)  using 
polynomial  approximation  of  K  given  in  Abramovitz 
and  Stegun  (1964) 

(10)  Determine  E,  A,  B,  f„  and  then  K  from  equation 

(12)  M  s 

(11)  Go  back  to  step  (4)  until  Kg  remains  constant 

(12)  Compare  the  obtained  K.  value  with  value  computed 
from  equation  (26).  If  different,  go  back  to  step 
(3) 

3 .  Results . 

The  wave  breaking  angle,  o^,,  as  computed  with  linear,  cnoidal- 
linear,  and  cnoidal  third-order  Stokes  approaches,  is  compared  with  the 
experimental  data  of  Kamphuis  (1963)  as  shown  in  Le  Mehaute  and  Koh  (1967) 
and  Figure  12.  These  are  apparently  the  only  data  on  a ^  and  are  ob¬ 
tained  for  a  single  bottom  slope,  S  =  0.1.  For  H0/L0  =  0.0175,  linear 
theory  gives  the  best  fit  to  the  data.  On  the  other  hand,  cnoidal 
linear  theory  provides  a  better  fit  for  the  larger  steepnesses,  H0/L0  = 
0.053  and  0.062.  Unfortunately,  it  is  not  possible  to  draw  any  de¬ 
finite  conclusions  from  the  limited  data.  However,  it  appears  that 
linear  theory  provides  the  best  estimate  of  wavelength  irrespective 
of  wave  steepness,  as  found  by  Eagleson  and  Dean  (Ippen,  1966).  Also, 
for  large  wave  steepness,  cnoidal  theory  predicts  the  wave  height  quite 
accurately  up  to  the  point  of  breaking. 

For  a  given  beach  slope  the  waves  break  at  increasing  relative 
depth  ratios,  d./L,  as  the  deepwater  steepness  increases.  For  large 
enough  H0/L0,  the  cnoidal  theory  is  no  longer  valid  since  it  predicts 
a  nonphysical  complex  wave  height  (Svendsen,  1974) .  The  critical 
deepwater  wave  steepness  for  which  the  cnoidal  theory  ceases  to  be 
valid  has  been  determined  for  five  different  bottom  slopes  (Table  2). 

This  critical  value  is  only  weakly  sensitive  to  the  magnitude  of  a0. 

For  H0/L0  greater  than  the  critical  value,  a  different  wave  theory 
such  as  Dean's  (1974)  or  third-order  Stokes  (Le  Mehaute  and  Koh,  1967) 
must  be  used. 

In  Figures  13  to  17,  the  variation  of  breaking  .,ave  angles 
versus  deepwater  wave  angle  with  Wq/Iq  and  S  as  parameters  is  de¬ 
picted  as  computed  using  the  onoidal-linear  theory.  Similar  results 
for  linear  theory  are  shown  in  Figure  18,  when  Wo/Lq  and  S  are  com¬ 
bined  into  a  single  parameter. 
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Figure  13.  Wave  breaking  angle,  aj,  as  a  fui 
deepwater  wave  angle,  a0,  and  dei 
steepness,  W0/L0,  on  a  bottom  slopi 
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Figure  14.  Wave  breaking  angle,  a^,  as  a  function  of 
deepwater  wave  angle,  a0,  and  deepwater 
steepness,  Ho/L^,  on  a  bottom  slope,  S  =  0: 


Figure  15.  Wave  breaking  angle,  a^,  as  a  function  of 
deepwater  wave  angle,  aQ,  and  deepwater 
steepness,  ti0/l0,  on  a  bottom  slope,  S  =  0 
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Figure  17.  Wave  breaking  angle,  a^,  as  a  tun* 
deepwater  wave  angle,  a0,  and  dee 
steepness,  \\0/L0,  on  a  bottom  slope 


Figure  18.  Prediction  of  wave  breaking  angles  using 
linear  theory. 


The  remaining  breaking  wave  characteristics  are  easily  determined 
as  shown  in  the  following  example: 


Given: 


Find : 
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(a)  Cnoidal -linear 


From  Figure  15  is  found 


% 


=  15°. 9 


Using  Snell 's  law 

h  sinab 

L  =  sina 
o  o 


0.426 


The  wave  height  can  be  found  from  the  breaking  criterion 

H.  H  -1/4  3/4 

jj-  -  0.76  cos  '  Oh  ( - ~)  Kr  =1.118 

no  L  S 

o 

Finally,  the  depth  of  breaking  is  computed  from  the  dispersion 
relationship,  which  is  solved  for 


or 
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(b)  Linear  theory 
H 

Determine  - -  0.148 

r?77 

0 

From  Figure  18 

=  12.2° 

Lb  sinab 


Snell's  law 


L  ~  sina 
o  o 


0.329 


Wave  height  fp  =  0.76  cos*^  ct^  ( - ^y) 


H,  -1/4  3/4 

K 


=  1.115 


Depth  at  breaking 


or 


k,  du  1  i„ 
b  b  =  y  In 


1  + 
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=  0.3417 


=  .0544  and  p-  =  p-  p"  =  0.0179 
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IV.  SUMMARY  AND  CONCLUSIONS 


It  appears  that  no  single  wave  theory  accurately  predicts  the  trans¬ 
formation  of  a  wave  from  deep  to  shallow  water.  In  shallow  water,  cnoidal 
theory  successfully  describes  the  shoaling  of  wave  height  but  overestimates 
the  wave  celerity  and  wavelength.  All  previous  studies  emphasize  an  accurate 
determination  of  wave  height,  but  little  attention  is  paid  to  the  wavelength. 

When  considering  longshore  currents  and  sediment  transport  the  angle 
of  wave  breaking  becomes  a  parameter  of  high  importance  and,  therefore,  also 
the  wave  refraction  process  which  is  intimately  connected  with  wavelength. 

In  this  study  a  "hybrid"  wave  theory  is  proposed,  consisting  of 
cnoidal  wave  height  and  linear  wavelength  transformation  in  the  regime  where 
cnoidal  theory  applies.  In  comparison  with  existing  data  this  theory  is 
found  to  predict  wave  height  and  wavelength  better  than  previous  theories, 
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especially  for  larger  wave  steepnesses.  More  definite  conclusions  must 
await  new  experimental  data.  This  basic  information  seems  to  be  missing 
in  the  existing  literature. 

Nomographs  are  presented  for  determining  the  breaking  wave 
characteristics  for  given  deepwater  characteristics  and  bottom  slope  for 
a  plane  beach. 

The  remaining  difficulties  encountered  in  the  shoaling  and  re¬ 
fraction  of  a  monochromatic  wave  are  further  amplified  when  dealing  with 
a  more  realistic  directional  spectrum  of  waves.  What  are  the  equivalent 
monochromatic  breaking  wave  characteristics  that  produce  the  same  sediment 
transport  or  longshore  current?  This  question  seems  to  be  one  of  the  more 
important  to  face  in  the  future  research. 
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APPENDIX 


Hyperbolic  Wave  Shoaling 


It  is  well  known  that  near  breaking  surface  gravity  water  waves  are 
highly  nonlinear  and  therefore  not  well  characterized  by  linear  theory. 

For  prediction  of  breaking  wave  characteristics  it  is  necessary  to  use 
higher  order  wave  theories  and,  in  particular,  a  higher  order  cnoidal  wave 
solution  could  appear  to  be  attractive  because  of  the  improved  fit  to  data 
found  when  using  first-order  cnoidal  theory.  Iwagaki  (1968)  derived  a 
practical  asymptotic  solution  to  second-order  cnoidal  waves  called 
"hyperbolic  waves."  This  approximation  takes  m  =  1  and  E(m)  =1  for 
K(m)  >  3  where  K(m)  and  E(m)  are  the  complete  elliptical  integrals  of  the 
first  and  second  kind. 

If  the  wave  height  transformation  is  written  in  the  usual  manner 
as  a  product  between  a  shoaling  and  refraction  coefficient 


H  =  H  K  .  K  (A-l) 

o  s  r 


(A-2) 


where 

d  =  water  depth 

L  =  (i  +  x2)  ,  deepwater  wavelength  (third-order  Stokes 
0  2ir  0  wave) 

=  unrefracted  deepwater  wave  height 

K  =  complete  elliptic  integral 

H  =  wave  height 
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dt  =  dfl  “  k  tl  +  II I  ^  1  =  water  dePth  under  trough 


The  parameter  \Q  is  determined  from 
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and  Iwagaki  found  by  empirical  curve  fitting 
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a  =  1.3,  n  =  2  and  m  *  j  for  j  <  0.55 

a  =  0.54,  n  =  \  and  m  =  1  for  y  >  0.55 
i  a 

Such  that  the  elliptic  integral,  K,  could  be  approximated  by 
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Finally,  in  hyperbolic  waves  the  wavelength  is  given  by 
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T~  =  /2it  (r~)  U  -  7-  (r-)  >  'T7T) 
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7  K  d 


n  ,H>nt"  rt  1  H,  5  II  ’ 

x  (1  -  a  (j)  }  {1  -  (1  +  k  d^  8  d 


(A-5) 


Again,  it  should  be  pointed  out  that  equations  (A-2)  to  (A-5)  are  only  valid 
for  K  >  3. 


To  evaluate  the  hyperbolic  wave  theory  the  shoaling  factor  Ks  is  com¬ 
pared  with  data  and  first-order  cnoidal  theory  (Brink-Kjaer  and  Jonsson, 
(1973).  The  results  are  shown  in  Figure  A-l. 


It  is  evident  that  the  results  given  by  the  hyperbolic  theory  are 
worse  than  those  obtained  using  linear  cnoidal  theory.  For  large  wave 
steepness  the  hyperbolic  theory  exhibits  a  decrease  in  Ks  near  breaking 
similar  to  fifth-order  Stokes  waves  which  is  also  due  to  nonhomogeneous 
convergence  of  the  perturbation  series.  Finally,  the  wavelength  predicted 
by  hyperbolic  theory  (not  shown  here)  is  quite  close  to  the  cnoidal  wave¬ 
length  which  is  in  itself  questionable  as  examined  earlier.  Further 
attempts  at  using  hyperbolic  wave  theory  to  predict  breaking  wave  character¬ 
istics  were  abandoned  because  of  the  above  shortcomings. 
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